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Abstract: The reaction mechanism of the gas-phase Pt atom with C3H8 has been 
systematically investigated on the singlet and triplet potential energy surfaces at 
CCSD(T)//BPW91/6-311++G(d, p), Lanl2dz level. Pt atom prefers the attack of primary 
over secondary C-H bonds in propane. For the Pt + C3H8 reaction, the major and minor 
reaction channels lead to PtC3H 6 + H 2 and PtCH 2 + C2H6, respectively, whereas the 
possibility to form products PtC2H4 + CH4 is so small that it can be neglected. The minimal 
energy reaction pathway for the formation of PtC3H 6 + H 2 , involving one spin inversion, 
prefers to start at the triplet state and afterward proceed along the singlet state. The optimal 
C-C bond cleavages are assigned to C-H bond activation as the first step, followed by 
cleavage of a C-C bond. The C-H insertion intermediates are kinetically favored over the 
C-C insertion intermediates. From C-C to C-H oxidative insertion, the lowering of 
activation barrier is mainly caused by the more stabilizing transition state interaction AE?mt 9 
which is the actual interaction energy between the deformed reactants in the transition state. 
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1. Introduction 

In recent years, the dehydrogenation of lower alkanes has gained great importance in natural and 
petroleum gas utilization [1,2]. Propane is a cheap and easily available raw material as it is produced 
through a number of petrochemical processes, while the propylene market demand is rapidly 
increasing. Accordingly, the dehydrogenation of propane is an interesting alternative route to 
propylene production. In the catalytic cracking of alkanes and catalytic dehydrogenation of alkanes, 
platinum based catalysts have received much attention. Elucidation of the role of isolated platinum 
units in heterogeneous catalysis can be aided by gas-phase study, which can provide insight into the 
intrinsic properties and reactivities of discrete and well-characterized catalytic species [3]. 

In particular, the activation of C-H and C-C bonds of propane by transition metals (neutral, cationic, 
or clusters) in the gas phase has been an active area of research that provides fundamental information 
on catalytic reaction mechanisms, kinetics, and thermodynamics [4,5]. Of the first-row transition metal 
series, for the activation of propane, the early members (Sc + [6], Ti + [7], and V + [8]) exhibit efficiency 
for the dehydrogenation of propane. Co + cation favors H 2 over CH 4 [9,10], whereas Fe + and Ni + 
cations favor CH 4 over H 2 [9,11]. Cr + cation does not show any efficiency for the activation of 
propane [12]. Of the second-row transition metal series, for Nb + [13], Mo + [12,14,15], and Rh + 
cations [16], the dehydrogenation of propane is efficient and the dominant process at low energies, 
whereas products resulting from both C-H and C-C cleavage processes are observable at high energies. 
Rh atoms are also effective for the H2 elimination from ethane and larger alkanes under kinetics 
technique experiment [17]. For the reaction of Ag + with propane, the dehydrogenation and formation 
of AgH + + R products are not observed, whereas the C-C bond cleavage is the predominant 
process [18]. Of the actinide ions, for the activation of propane, Th + , Pa + and U + cations are efficient 
for the dehydrogenation of propane. For the other cations (Np + , Pu + , Am + , and Cm + ), no reactions are 
observed experimentally [19]. For Th and U , both C-H and C-C cleavage products are effectively 
observed [20]. Finally, for [(MgO)„] + clusters, the higher reactivity with propane is not specific to 
[(MgO) 2 ] + , but has been also observed for (MgO) + [21,22]. 

Concerning the 5d-series transition metal Pt (neutral, cationic, or clusters), the reactions with linear 
alkanes have been extensively explored by means of diverse experimental and theoretical 
methods [23-29]. It is reported that the transition metal Pt (neutral, cationic, or clusters) are the 
efficient C-H insertion agents [23-29]. However, as far as we know, although a few investigations of 
C-H insertion processes have focused on the neutral Pt atom [23-27], the C-C insertion has hardly 
been investigated in the gas phase. 

In the present study, a complete mechanism of neutral Pt atom with propane along with both the 
C-H and C-C bond activation processes is investigated, which is necessary to enable us to determine 
the crucial steps and to either block or enhance particular steps to steer the reaction in the desired 
direction. The goals of the present investigation are as follows: (a) to provide reliable structures and 
chemically accurate energetics of the reactants, intermediates, transition states (TSs), and products; 
(b) to elucidate the rate-determining step and the selectivity-controlling step; and (c) to gain a 
better understanding of the preference of reaction pathway. Particularly, to shed some light on the 
intrinsic reactivity of platinum atom toward the propane activation, the trends in reactivity and 
competition among the C-H and C-C bond cleavage mechanisms are analyzed using the activation 



Int. J. Mol Sci. 2012, 13 



9280 



strain model of chemical reactivity [30,31]. The potential energy profiles on the singlet and triplet 
states are investigated, because spin crossing is often involved in the transition metal-containing 
reactions [7,15,20]. 

2. Computational Details 

All calculations were carried out with the Gaussian 03 program [32]. Full geometry optimizations 
were run to locate all of the stationary points and TSs on the singlet and triplet potential energy 
surfaces (PESs) for the reactions of Pt atom with C3H8, using the BPW91 [33,34] method with 
6-31 1++G(d, p) basis set for the carbon and hydrogen atoms [35,36], and the Lanl2dz basis set and the 
corresponding effective core potential (ECP) for platinum [37], namely BPW91/6-311++G(d, p), 
Lanl2dz. Meantime, the stability of the density function theory (DFT) wavefunction was tested [38,39]. 
If an instability was found, the wavefunction was reoptimized with appropriate reduction in 
constraints, and the stability tests and reoptimizations were repeated until a stable wavefunction was 
found [38,39]. Harmonic frequency calculations were run to characterize stationary points and to take 
corrections of zero-point energy (ZPE) into account. The intrinsic reaction coordinate (IRC) method 
was performed to track minimum energy paths from transition structure to the corresponding 
minima [40,41]. The dominant occupancies of natural bond orbitals for some species have been 
analyzed with the help of the natural bond orbital (NBO) analysis [42,43]. To further determine 
electron correction energies, the single-point calculation of various species based on the optimized 
BPW91/6-311++G(d, p), Lanl2dz geometries were then refined using CCSD(T) [44] method 
with the same basis sets, namely CCSD(T)//BPW91/6-311++G(d, p), Lanl2dz. Unless 
otherwise mentioned, all energies are relative to the ground-state reactants [Pt( 3 D) + C3H8] at the 
CCSD(T)//BPW91/6-311++G(d, p), Lanl2dz level, including ZPE correction obtained at the 
BPW91/6-31 1++G(d, p), Lanl2dz level. 

3. Results and Discussion 

Considering the present system, the BPW9 1/6-31 lG++(d, p), Lanl2dz level is suitable to reproduce 
experimental values of geometrical parameters of H 2 , Pt 2 , and PtH diatomic molecules [45]. 
Furthermore, the single-point calculation of various species were then refined at the 
CCSD(T)//BPW91/6-311++G(d, p), Lanl2dz level. Thereupon, the present theoretical method of 
CCSD(T)//BPW91/6-3 1 1++G(d, p), Lanl2dz should be appropriate and reliable for the Pt + C 3 H 8 systems. 

In this work, we will mainly discuss the following reactions of Pt atom with C3H8: 

Pt + C 3 H 8 ^PtC 3 H 6 + H 2 (1) 

Pt + C 3 H 8 ^PtCH 2 + C 2 H 6 (2) 

Pt + C 3 H 8 ^PtC 2 H 4 + CH 4 (3) 

The above three product channels were divided into five sections: (i) C-H bond activation: 
dehydrogenation, (ii) C-C bond activation: deethanization, (iii) C-C bond activation: demethanation, 
(iv) comparison of C-H with C-C bond activation, and (v) activation strain analysis of the direct C-H 
and C-C bond cleavage. 
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3.1. C-H Bond Activation: Dehydrogenation 

For the dehydrogenation of C3H8 by Pt atom, the reaction pathway and the optimized geometric 
structures of various species are depicted in Scheme 1 . 

The triplet state Pt atom, 3 D(dV), is the ground state. The singlet state Pt atom, ^(d 10 ), lies 
59.4 kJ-mor 1 above the ground triplet state 3 D(dV), in good agreement with the estimated value of 

1 13 

56.6 kJ-mol [26]. The superscript prefixes " " and " " will be used to indicate the singlet and triplet 

states, respectively. As depicted in Scheme 1, there are two primary (tans-C-H^ and cis- C-H (1) in the 

trans and cis position with respect to the CH3 group, respectively) and one secondary (C-H (2) ) C-H 

bonds in propane. Then, with regard to the initial interaction between Pt atom and C3H8, three 

molecular complexes are considered: (i) Pt atom attacking the H-end of primary tans-C-H^ 

(l-PtCsHg), (ii) Pt atom attacking the H-end of secondary C-H (2) approaching to Pt atom (2-PtC3H 8 ), 

(iii) Pt atom attacking the H-end of primary cis-C-H^ (3-PtC3H 8 ). 

As discussed earlier, Pt atom has a triplet ground state ( D) with excitation energy of 59.4 kJ-mol 

to the lowest singlet state ( 1 S). Considering the initial interaction of Pt atom with C3H8, only the triplet 

ground state 3 l-PtC 3 H 8 , 3 2-PtC3H 8 , and 3 3-PtC3H 8 molecular complexes are obtained, whereas we 

failed to locate the corresponding ones on the singlet PES despite extensive attempts. For 3 l-PtC3H 8 , 

3 2-PtC 3 H 8 , and 3 3-PtC 3 H 8 , the BSSEs [46] by BPW91 are 10.1, 10.5, and 10.3 kJ-mol" 1 , and the 

complexation energies corrected by BSSEs are calculated to be 1.6, 15.7, and 3.8 kJ-mol 1 relative to 

the reactants Pt( D) + C3H 8 , respectively. It is shown that the complex stability increases along 

3 l-PtC 3 H 8 < 3 3-PtC 3 H 8 < 3 2-PtC 3 H 8 . For 3 l-PtC 3 H 8 , 3 2-PtC 3 H 8 , and 3 3-PtC 3 H 8 , the C-H bond close to 

Pt atom is elongated to 1.147, 1.164 and 1.147 A from the 1.100, 1.103, and 1.101 A of free C 3 H 8 , 

while there is a short Pt-H distance of 2.016, 1.955, 2.023 A, respectively, indicating some molecular 

interaction between Pt atom and C3H 8 . The minimal energy reaction pathway (MERP) may start at the 

1.1. 1. 
triplet molecular complexes ( l-PtCyig, 2-PtC3H 8 , and 3-PtC3H 8 ) from the corresponding ground 

triplet reactants. 

As shown in Scheme 1, from these molecular complexes (l-PtCsHg, and 2-PtC3H 8 , 3-PtC3H 8 ), 
the C-H bond cleavage may lead to the dehydrogenation product PtC3H 6 (Pt-propene) + H 2 and 
Pt(CH2)3 (Pt-cyclopropane) + H2. The change of Gibbs free energies (AG29 8 ) for the reactions of 
Pt( 3 D) + C 3 H 8 ^ 1 PtC 3 H 6 (Pt-propene) + H 2 and 1 Pt(CH 2 ) 3 (Pt-cyclopropane) + H 2 are calculated to 
be -111.6 and -104.5 kJ-mol 1 , respectively. Thereby, the dehydrogenation of C3H 8 is 
thermodynamically favorable. 

For the formation of PtC3H 6 (Pt-propene) + H 2 , there are five reaction pathways beginning at the 
three kinds of molecular complexes (two from l-PtC3H 8 , two from 2-PtC3H 8 , and one from 3-PtC3H 8 ), 
respectively, as shown in Scheme 1. Alternatively, for the formation of Pt(CH 2 )3 (Pt-cyclopropane) + H 2 , 
there is an unique reaction pathway starting from 3-PtC3H 8 . 
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Scheme 1. The reaction pathway and the optimized geometric structures of various species in the dehydrogenation of C3H8 by Pt atom, 
through (a) l-PtCsHg; (b) 2-PtC3Hs, and (c) 3-PtC3Hs. Bond lengths are reported in A and bonds angles in degree. Relative energies (kJ-moi ) 
for the corresponding species relative to Pt( 3 D) + C 3 H 8 at the CCSD(T)//BPW9 1/6-31 1++G(d, p), Lanl2dz level are shown. 




(b) 
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First, from l-PtCsHg, the initial primary C-H bond oxidative insertion via TS1 or TS2 leads to the 
a-complex intermediate l-HPtC3H 7 or c/s-l-HPtC3H 7 , respectively. From l-HPtC3H 7 , a 

1 .2- dehydrogenation process takes place via four-center transition state TS3, resulting in a dihydrogen 
propene complex (H2)PtC3H6. Finally, the molecular complex (H2)PtC3H6 reductively eliminates H 2 , 
leaving PtC3H 6 behind. Alternatively, from cis-l-HP^Hy, a 1 ,2-dehydrogenation process occurs via 
five-center transition state TS4 directly, leading to the dissociation products PtC3H 6 + H 2 . 

Second, from 2-PtC3H 8 , the initial secondary C-H bond oxidative insertion via TS5 yields the 
a-complex intermediate HPtCH(CH 3 ) 2 . From HPtCH(CH 3 ) 2 , there are two reaction pathways for the 
formation of PtC3H 6 + H 2 . On the one hand, a 1 ,2-dehydrogenation process takes place via four-center 
transition state TS6, producing the dihydrogen propene complex (H 2 )PtC3H 6 . On the other hand, from 
HPtCH(CH 3 ) 2 , a 1,2-dehydrogenation process occurs via five-center transition state TS7, directly 
resulting in the dissociation products PtC3H 6 + H 2 . 

Third, from 3-PtC3H 8 , the initial primary C-H bond oxidative insertion via TS8 or TS9 generates a 
a-complex intermediate 3-HPtC3H 7 or c/s-3-HPtC3H 7 , respectively. From 3-HPtC3H 7 , a 

1.3- dehydrogenation process takes place via five-member transition state TS10, generating the 
dihydrogen metallacycle molecular complex H 2 Pt(CH 2 )3. Last, the molecular complex dissociates into 
Pt(CH 2 ) 3 + H 2 . The structure of metallacycle Pt(CH 2 ) 3 is similar to that of Sc(CH 2 ) 3 + [6], TiC 3 H 6 + [7], 
and NiC4H 8 + [47]. Alternatively, from c/s-3-HPtC3H 7 , a 1,2-dehydrogenation process occurs via 
five-center transition state TS11, producing the dihydrogen propene complex (H 2 )PtC3H6. Finally, 
(H 2 )PtC3H 6 reductively eliminates H 2 , leaving PtC3H 6 behind. 

For the formation of PtC3H 6 + H 2 , from l-PC^Hg, the MERP should proceed via the minimal 
energy crossing point (MECP) between 1 l-HPtC3H 7 and 3 l-HPtC3H 7 , with the energy height of the 
highest point (EHHP) of 61.2 kJ-mol" 1 at 3 TS1. From 2-PtC 3 H 8 , the MERP should proceed via MECP 
between 1 HPtCH(CH 3 ) 2 and 3 HPtCH(CH 3 ) 2 , with the EHHP of 69.5 kJ-mol" 1 at 3 TS5. From 3-PtC 3 H 8 , 
the MERP should progress via MECP between 1 c^-3-HPtC 3 H 7 and 3 c/s-3-HPtC 3 H 7 , with the EHHP of 
162.1 kJ-mol" 1 at TS1 1. Since the EHHP at 3 TS1 from l-PtC 3 H 8 is the lowest among the three reaction 
channels, this reaction channel for the formation of PtC3H 6 + H 2 is the most feasible kinetically. 
Furthermore, these results reveal a high preference of Pt atom for the attack of primary C-H bonds in 
propane, which is analogous to that of MgO + cation for the attack of alkanes [21]. This feature 
represents a notable distinction of the transition-metal atom from various transition-metal oxide 
cations, which show a clear preference for the attack of secondary C-H bonds [21]. 

For the formation of Pt(CH 2 ) 3 + H 2 , from 3-PtC 3 H 8 , the MERP should advance via MECP between 
^-HPtCsHv and 3 3-HPtC 3 H 7 , with the EHHP of 65.5 kJ-mol" 1 at 3 TS8. 

Moreover, the 1 l-HPtC 3 H 7 , Wl-HPtC 3 H 7 , 1 HPtCH(CH 3 ) 2 , 1 3-HPtC 3 H 7 , W3-HPtC 3 H 7 , and 
1 (H 2 )PtC 3 H 6 intermediates lie -182.5, -181.6, -192.3, -180.4, -183.3, and -202.9 kJ-mol" 1 in a deep 
energetic well on each MERP, respectively. Then, these intermediates are thermodynamically favored 
in the dehydrogenation of C3H 8 . For the intermediates containing -PtH and -Pt-alkyl moieties 
(H-HPtCsHv, Wl-HPtC 3 H 7 , 1 HPtCH(CH 3 ) 2 , ^-HPtCsHy, and W3-HPtC 3 H 7 ), the NBO results 
show that a complete a-bond has been formed both in Pt-H and in Pt-C. 

From Pt + C 3 H 8 to the C-H insertion intermediates ( 1 l-HPtC 3 H 7 , Wl-HPtC 3 H 7 , 1 HPtCH(CH 3 ) 2 , 
1 3-HPtC3H 7 , and 1 c/.s"-3-HPtC3H 7 ), only the triplet molecular complexes and the triplet TSs are 
obtained, while we failed to gain the corresponding singlet ones, despite extensive attempts. 
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Furthermore, for Pt atom (10 valence electrons), Pt^S) singlet state, has an empty orbital and five 
doubly occupied nonbonding orbitals, whereas Pt( 3 D) triplet state has all of its s and d valence orbitals 
occupied, with four doubly occupied nonbonding orbitals, and two singly occupied nonbonding 
orbitals. That is to say, the bonding capacity of Pt^S) singlet state to C3H8 is stronger than that of 
Pt( D) triplet state. Then, the binding of Pt to C3H8 in the C-H insertion intermediates (l-HPtC3H 7 , 
c/s-l-HPtC 3 H 7 , HPtCH(CH 3 ) 2 , 3-HPtC 3 H 7 , and c/s-3-HPtC 3 H 7 ) inverts the energies of the singlet and 
triplet states from the Pt atom. Therefore, the ground state of the C-H insertion intermediates is the 
singlet state, as depicted in Scheme 1. The reaction goes forward from the excited state reactants Pt^S) 
+ C 3 H 8 to the C-H inserted intermediates, without energy barrier. This can be ascribed to the fact that 
the singlet state Pt( 1 S) has an empty orbital, which should greatly facilitate the interaction with 
propane and its bond activation, leading smoothly to the formation of two covalent bonds to H and 
C 3 H 7 . This phenomenon has also appeared in the analogous Pt + CH 4 system [26]. 

3.2. C-C Bond Activation: Deethanization 

For the deethanization of C 3 H 8 by Pt atom, the reaction pathway and the optimized geometric 
structures of various species are depicted in Scheme 2. The change of Gibbs free energies (AG 2 9 8 ) for 
the reaction of Pt( 3 D) + C 3 H 8 ^ 1 PtCH 2 + C 2 H 6 is calculated to be -45.9 kJ-mol" 1 . Thereupon, the 
deethanization of C 3 H 8 is thermodynamically favorable. Then, it is necessary to discuss kinetically the 
above reaction infra. 

As shown in Scheme 2, for the deethanization of C 3 H 8 by Pt atom, there are seven reaction 
pathways, three from l-PtC 3 H 8 , one from 2-PtC 3 H 8 , and three from 3-PtC 3 H 8 . These seven reaction 
pathways are separated into two kinds of reaction pathways, one through the initial C-C bond direct 
cleavage, and another through the a-complex assisted C-C a-bond metathesis. For simplicity, we will 
primarily discuss the three reaction pathways from l-PtC 3 H 8 infra, which are analogous to those from 
2-PtC 3 H 8 and 3-PtC 3 H 8 . The nuances in energy mainly stem from their configuration differences 
among the three kinds of reaction channels. 

As mentioned earlier, from l-PtC 3 H 8 , there are three reaction pathways for the PtCH 2 + C 2 H 6 
formation. For the first reaction pathway, Pt atom firstly inserts the C-H bond via five-member TS2, 
resulting in the intermediate c/s-l-HPtC 3 H 7 . Then, from cis-l-HPtC 3 H 7 , a a-complex assisted C-C 
a-bond metathesis takes place via a four-member TS12 with both 1,3-H migration and C-C cleavage, 
yielding the molecular complex, C 2 H 6 PtCH 2 . Finally, C 2 H 6 PtCH 2 releases C 2 H 6 molecule, leaving 
PtCH 2 behind. The MERP should proceed via the MECP between Wl-HPtC 3 H 7 and Wl-HPtC 3 H 7 , 
with the highest energy requirement (HER) of 273.2 kJ-mol 1 at the 1 c/.s"-l-HPtC 3 H 7 ^ 1 TS12 reaction 
step and the EHHP of 91.6 kJ-mol" 1 at l TS\2. 

For the second and third reaction pathways from l-PtC 3 H 8 , there are two reaction pathways to 
produce c/s-CH 3 PtC 2 H 5 . That is, Pt atom directly inserts the C-C bond via three-member TS13, 
leading to the intermediate c/s-CH 3 PtC 2 H 5 . Alternatively, Pt atom firstly inserts the C-H bond via the 
three-member TS1, resulting in the intermediate l-HPtC 3 H 7 . Next, from l-HPtC 3 H 7 , a a-complex 
assisted C-C a-bond metathesis occurs via a four-member TS14 with both 1,2-H shift and C-C bond 
cleavage, also yielding c/s-CH 3 PtC 2 H 5 . Then, from CH 3 PtC 2 H 5 , the a-complex assisted C-H a-bond 
metathesis takes place via a four-member TS15 with 1,3-H shift, yielding the molecular complex 
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C2H6PtCH2. As mentioned earlier, C2H6PtCH2 releases C2H6 molecule, staying PtCH2 behind. For the 
two reaction pathways, each MERP should advance via the MECP between 1 c^-CH 3 PtC2H 5 and 
3 c/s-CH 3 PtC 2 H 5 , with the HER of 299.0 kJmor 1 at the 1 c^-CH 3 PtC 2 H 5 ^ 1 TS15 reaction step. The 
two MERPs involve the EHHP of 156.1 and 120.3 kJ-mol" 1 at 3 TS13 and ^SIS, respectively. 

Scheme 2. The reaction pathway and the optimized geometric structures of various species 
in the deethanization of C3H8 by Pt atom, through (a) l-P^Hg and 2-PtC3H 8 , and 
(b) 3-PtC3Hg. Bond lengths are reported in A and bonds angles in degree. Relative energies 
(kJ-mol 1 ) for the corresponding species relative to Pt( 3 D) + C3H8 at the 
CCSD(T)//BPW91/6-31 1++G(d, p), Lanl2dz level are shown. 




3-HPtC 3 H 7 ! TS17 



(b) 
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Comparing these three reaction pathways, one can conclude that the reaction pathway via TS2 
involving the first C-H cleavage and via TS12 involving C-C cleavage with synchronous 1,3-H 
migration is the gross MERP for the l-PC^Hg— >PtCH 2 + C2H6 reaction, owing to its comparatively 
low HER (273.2 vs. 299.0 kJ-mol" 1 ) and low EHHP (91.6 vs. 120.3 and 156.1 kJ-mol" 1 ), with the 
rate-determining step of 1 c^-l-HPtC 3 H7^ 1 TS12^ 1 C2H 6 PtCH2. 

From 2-PtC3H 8 for the formation of PtCH 2 + C2H6, only one reaction pathway is obtained, which 
includes the initial C-C bond cleavage and a-complex assisted C-H a-bond metathesis. This reaction 
pathway is homologous to that via the initial C-C cleavage from l-PtC3H 8 with the HER of 

299.0 kJ-mol" 1 and the EHHP of 156.1 kJ-mol" 1 at 3 TS13. 

From 3-PtC3Hg for the formation of PtCH2 + C2H6, there are also three reaction pathways. These 

three reaction pathways are similar to those from l-PtCsHg. The reaction pathway of a-complex 

3 1 

assisted C-C a-bond metathesis via TS9 and TS19 is kinetically most preferable in the three reaction 
pathways, because of its lowest HER (274.9 vs. 298.7 kJ-mol" 1 ) and lowest EHHP (91.6 vs. 120.3 and 

156.1 kJ-mol" 1 ). 

In summary, for the formation of the C-C cleavage products PtCH 2 + C2H6, the optimal pathway 
proceeds through the a-complex c/s-l-HPtC3H 7 or c/s-3-HPtC3H 7 from initial C-H bond cleavage, 
which assists the C-C a-bond metathesis. This reactivity mode is also complementary for the classical 
reactivity picture through the direct C-C cleavage intermediate (M = Fe + [48] and Ta + [49]). 

3.3. C-C Bond Activation: Demethanation 

For the demethanation of C3H8 by Pt atom, the reaction pathway and the optimized geometric 
structures of various species are depicted in Scheme 3. The change of Gibbs free energies (AG298) for 
the reactions of Pt(3D) + C 3 H 8 ^lPtC 2 H 4 + CH 4 are calculated to be -143.6 kJ-mol" 1 , which is 
thermodynamically favorable. Afterwards, we will discuss the kinetics of the above reaction infra. As 
shown in Scheme 3, there are four kinds of reaction pathways, which are through c/s-CH 3 PtC2H 5 , 
CH 3 PtC 2 H 5 , HPtCH(CH 3 ) 2 , and C 2 H 4 PtH(CH 3 ) intermediates, respectively. 

First, from c/s-CH 3 PtC2H 5 , a a-complex assisted C-H a-bond metathesis takes place via a 
four-member TS20 with 1,4-H shift, yielding the molecular complex CH 4 PtC2H 4 . The molecular 
complex CH 4 PtC2H 4 releases CH 4 molecule, leaving PtC2H 4 behind. Through c/s-CH 3 PtC2H 5 , the 
MERP should go forward via the MECP between 1 l-HPtC 3 H 7 and 3 l-HPtC 3 H 7 , with the HER of 
319.0 kJ-mol" 1 at the 1 ds-CH 3 PtC 2 H 5 ^ 1 TS20 reaction step and the EHHP of 140.3 kJ-mol" 1 at 1 TS20. 

Second, from CH3PtC2H 5 , a a-complex assisted C-H a-bond metathesis occurs via a four-member 
TS27 with 1,3-H shift, yielding a molecular complex CHCH 3 PtCH 4 . Then, the molecular complex 
CH 4 PtCHCH 3 sets a CH 4 molecule free, leaving PtCHCH 3 behind. Next, from PtCHCH 3 , 1,2 H shift 
occurs via a four-member TS22, staying PtC2H 4 behind. Through CH3PtC2H 5 , the MERP should go 
forward via the MECP between 1 3-HPtC 3 H 7 and 3 3-HPtC 3 H 7 , with the HER of 225.8 kJ-mol" 1 at the 
1 3-HPtC 3 H 7 ^ 1 TS17 reaction step and the EHHP of 68.9 kJ-mol" 1 at 1 TS22. 

Third, from HPtCH(CH3)2, a a-complex assisted C-C a-bond metathesis occurs via a four-member 
TS21 with 1,3-H shift, also leading to the molecular complex CHCH 3 PtCH 4 . The MERP should go 
forward via the MECP between 1 2-HPtC 3 H 7 and 3 2-HPtC 3 H 7 , with the HER of 277.4 kJ-mol" 1 at the 
1 HPtCH(CH 3 ) 2 -> 1 TS21 reaction step and the EHHP of 85.1 kJ-mol" 1 at 1 TS21. 
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Scheme 3. The reaction pathway and the optimized geometric structures of various species in the demethanation of CsH 8 by Pt atom, through 
(a) l-PtC3Hs and 2-PtC3Hs, and (b) 3-PtC3H8. Bond lengths are reported in A and bonds angles in degree. Relative energies (kJ-mol - ) for the 
corresponding species relative to Pt( 3 D) + C 2 H 6 at the CCSD(T)//BPW9 1/6-31 1++G(d, p), Lanl2dz level are shown. 
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Scheme 3. Cont. 
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Fourth, from 3-HPtC3H 7 , an oxidative insertion of C-C bond to the platinum center takes place via a 
four-member TS23, producing a methyl hydrid complex C2H 4 PtH(CH 3 ). Then, from C2H 4 PtH(CH 3 ), 
1,2-H shift occurs, yielding the molecular complex CH 4 PtC2H 4 . Last, the molecular complex 
CH 4 PtC2H 4 sets a CH 4 molecule free, leaving PtC2H 4 behind. The MERP should go forward via the 
MECP between ^-HPtCsHy and 3 3-HPtC 3 H 7 , with the HER of 143.1 kJ-mol" 1 at the 
1 3-HPtC 3 H 7 ^ 1 TS23 reaction step and the EHHP of 65.5 kJ-mol" 1 at 3 TS8. 

Comparing these four kinds of reaction pathways, one can see that the reaction pathway starting at 
the 3-PtC3H 8 involving the crucial intermediate C2H 4 PtH(CH 3 ) is the most optimal MERP for the 
Pt + C 3 H 8 ^PtC 2 H 4 + CH 4 reaction, thanks to the lowest HER (143.1 vs. 319.0, 225.8, and 
277.4 kJ-mol" 1 ) and lowest EHHP (65.5 vs. 140.3, 85.1, and 68.9 kJ-mol" 1 ), with the rate-determining 
step of 1 3-HPtC 3 H 7 ^ 1 TS23^C 2 H 4 PtH(CH 3 ). Thereby, the optimal pathway proceeds through the 
a-complex 3-HPtC 3 H 7 from initial C-H bond cleavage, which assists the C-C a-bond metathesis. That 
is to say, the optimal C-C bond cleavages are assigned to C-H bond activation as the first step, 
followed by cleavage of a C-C bond. This reactivity mode is complementary for the classical reactivity 
picture through the direct C-C cleavage intermediate (M = Fe + [48], and Ta + [49]). 

3.4. Comparison of C-H with C-C Bond Activation 

As shown in Schemes 1-3, the C-H insertion intermediates ( 1 l-HPtC 3 H 7 , 1 c/.s"-l-HPtC 3 H 7 , 
1 HPtCH(CH 3 ) 2 , 1 3-HPtC 3 H 7 , W3-HPtC 3 H 7 ) and the C-C insertion intermediates ^CHsPtC^s, 
1 c^-CH 3 PtC 2 H 5 , 1 C 2 H 4 PtH(CH 3 ), and 1 CH 4 PtC 2 H 4 ) deposit in a deep well, respectively. It is indicated 
that these intermediates are thermodynamically preferred. For the formation of the C-H and the C-C 
insertion intermediates, the corresponding MERP should involve the HER of about 60-70 and 
140-230 kJ-mol 1 , respectively. Thereby, the C-H insertion intermediates are kinetically favored, 
while the C-C insertion intermediates are kinetically hindered by energy barriers. These results are in 
qualitative agreement with the experimental results, in which the C-H insertion product is 
experimentally observed and the C-C insertion product is not formed in observable quantity in 
Pt + C 2 H 6 system [27]. 

For the formation of C-C bond cleavage intermediates 1 CH 3 PtC2H 5 and 1 c/^-CH 3 PtC2H 5 , one can 
see that the reaction pathways of the direct C-C activation via 3 TS16and 3 TS13 are inferior to those of 
the a-complex assisted C-C a-bond metathesis via 1 TS17 and X TSU from 3-HPtC 3 H 7 and l-HPtC 3 H 7 , 
respectively, because of their higher EHHP (156.1 vs. 65.5 and 61.2 kJ-mol 1 ). This is reminiscent of 
the important role of a-complex assistance for the C-C a-bond metathesis. In other words, the direct 
C-C bond activation is associated with a sizable barrier, which would prohibit this channel. 

A glance to the reaction pathways shown in Schemes 1-3 reveals that two kinds of a-complexes 
(l-HPtC 3 H 7 , HPtCH(CH 3 ) 2 , and 3-HPtC 3 H 7 ) and (cfe-l-HPtC 3 H 7 and cw-3-HPtC 3 H 7 ) from initial C-H 
bond cleavage are crucial for the selective formation of the final C-H and C-C cleavage products. 

First, from the identical intermediate l-HPtC 3 H 7 , the reaction step of 1 l-HPtC 3 H 7 ^ 1 TS3^ 1 (H 2 )PtC 3 H 6 
is competitive with that of H-HPt^H^TSM^ WCH 3 PtC 2 H 5 . Because *TS3 lies 160.9 kJ-mol" 1 
below ^SH, (H 2 )PtC 3 H 6 is selectively preferred, whereas c/s-CH 3 PtC2H 5 is selectively hampered. In 
other words, from l-HPtC 3 H 7 , the dehydrogenation process dominates. 
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Second, from the identical intermediate c/s-l-HPtC 3 H 7 , the reaction step of 
1 c^-l-HPtC 3 H7^ 1 TS12^ 1 C 2 H 6 PtCH2 is competitive with that of 1 c^-l-HPtC3H7^ 1 TS4^PtC 3 H 6 + H 2 . 
Since *TS12 locates 57.3 kJ-mol 1 below 1 TS4, C2H 6 PtCH 2 is selectively favored. That is to say, from 
c/s-l-HPtC 3 H 7 , the deethanization process predominates. 

Third, from the identical intermediate HPtCH(CH 3 ) 2 , these reaction steps of HPtCH(CH 3 ) 2 ^ 1 TS6^ 
(H 2 )PtC 3 H 6 , HPtCH(CH 3 ) 2 ^ 1 TS7^PtC 3 H 6 + H 2 , and HPtCH(CH 3 ) 2 ^ 1 TS21^CH 3 CHPtCH 4 are 
competitive. As l TS6 lies 276.5 and 210.7 kJ-mol" 1 below l TS7 and 1 TS21, respectively, (H 2 )PtC 3 H 6 is 
selectively preferred. Then, from HPtCH(CH 3 ) 2 , the dehydrogenation process dominates. 

Fourth, from the identical intermediate 3-HPtC 3 H 7 , these reaction steps of 3-HPtC 3 H 7 — ^TSIO— ► 
(H 2 )Pt(CH 2 ) 3 , 3-HPtC 3 H 7 ^ 1 TS17^CH 3 PtC 2 H 5 , and 3-HPtC 3 H 7 ^ 1 TS23^C 2 H 4 PtH(CH 3 ) are 
competitive. Because ^SIO locates 117.0 and 34.3 kJ-mol 1 below *TS17 and ^823, respectively, 
(H 2 )Pt(CH 2 ) 3 is selectively favored, whereas CH 3 PtC 2 H 5 and C 2 H 4 PtH(CH 3 ) are selectively hindered. 
Thereby, from 3-HPtC 3 H 7 , the dehydrogenation process dominates. 

Last, from the identical intermediate c/s-3-HPtC 3 H 7 , these reaction steps of c/s-3-HPtC 3 H 7 — ^TSl 1 
-^PtC 3 H 6 + H 2 and c^-3-HPtC 3 H 7 ^ 1 TS19^ 1 C 2 H 6 PtCH 2 are competitive. Because *TS19 lies 
70.5 kJ-mol" 1 below ^Sll, 1 C 2 H 6 PtCH 2 is selectively preferred. Therefore, from c/s-3-HPtC 3 H 7 , the 
deethanization process predominates. 

In summary, once the a-complex [l-HPtC 3 H 7 , or HPtCH(CH 3 ) 2 , or 3-HPtC 3 H 7 ] is formed, the 
major reaction channel results in the dehydrogenations products. Alternatively, as far as the a-complex 
[c/s-l-HPtC 3 H 7 or c/s-3-HPtC 3 H 7 ] is concerned, the major reaction channel leads to the deethanization 
products. Besides, the demethanation process is kinetically ruled out. 

To estimate quantitatively the reactivity and selectivity for the two kinds of products [PtC 3 H6 + H 2 
and PtCH 2 + C 2 H 6 ], the rate constants have been evaluated according to conventional transition state 
theory (TST) [50], including tunneling correction based on Winger's formulation [51]. The formation 
of rate constant k(T) including tunneling correction coefficient #"(T) in transition state theory is 

given by 

*(T) = * / (T)*at(T) (4) 
The rate constant k'{T) is simply given by 

k T — 

where k B is the Bolzmann constant, h is the Planck constant, T is thermodynamic temperature, c° is 
standard concentration, and AG^ is Gibbs free energy. The tunneling correction coefficient /c(T) is 
written in the form of 

, . 1 hv* 2 

(6) 

where k B is the Bolzmann constant, h is the Planck constant, T is thermodynamic temperature, and v* 
is the imaginary frequency of the unbound normal mode at the saddle point. The branching ratio («i) of 
product / is calculated by 



1 


hv* 


24 


k B T 
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_ *,(T) 
Z^.(T) 



a =^ (7) 



where ^ is the rate constant of product /. 

From the identical reactants Pt + C3H8, the formation of PtC3H6 + H2 and PtCH2 + C2H6 are 

o 11 

competitive, while their selectivity-controlling steps are Pt( D) + C3H8— » TS1— ► l-HPtC3H 7 and 
Pt( 3 D) + C 3 H 8 -> 3 TS9^ 3 c/s-3-HPtC 3 H 7 on their MERPs, respectively. Thereby, the rate constants 

o 11 

were taken into account, where Pt( D) + C3H8 were taken as reactants, while TS1 and TS9 served as 
TSs, respectively. The rate constants for the formation of 3 l-HPtC3H 7 (£1) and 3 c/s-3-HPtC3H 7 (£ 2 ) 
calculated over 300-1100 K temperature range can be fitted by the following expressions 
(in dm^mol ^s 1 ): 

k x = 6.65x10 s exp(-65,855/i?r) (8) 
k 2 = 2.85x10 s exp(-72, 999 1 RT) (9) 

The branching ratios for the formation of PtC3H6 + H2 and PtCH2 + C2H6 are calculated to be 
97.7-83.6% and 2.3-16.4%, respectively, over 300-1100 K temperature range. In other words, the 
dehydrogenation channel is predominant, and the deethanization channel is minor, while the 
demethanation channel is ruled out. 



3.5. Activation Strain Analysis of the direct C-H and C-C Bond Cleavage 

To gain insight into how the Pt atom affects the activation barriers of the initial C-H and C-C bond 
cleavage, i.e., insight into how this effect depends on the nature of concomitant geometrical 
deformation and electronic structures of Pt and C3H8, the trends in reactivity and competition among 
the initial C-H and C-C bond mechanisms are analyzed using the activation strain model of chemical 
reactivity [30,31]. In this model, activation energies AEf of the TS are divided into the activation strain 
Ai^strain and the stabilizing TS interaction AE^[ nt : AEf = AE^ stm [ n + AE^- mi . The activation strain AE^ stm [ n 
is the strain energy associated with deforming the reactants from their equilibrium geometry to the 
geometry they adopt in the TS. The TS interaction AE^{ nt is the actual interaction energy between the 
deformed reactants in the TS [30,31]. 

The results of the activation strain analysis are listed in Table 1. The activation energy AE* 
increases from the initial C-H cleavage TSs ( 3 TS1, 3 TS2, 3 TS5, 3 TS9, and 3 TS8) of 60-70 kJ-mol" 1 to 
the initial C-C cleavage TSs ( 3 TS13 and 3 TS16) of -160 kJ-mol -1 . The activation strain A£^ str ain 
decreases from the intial C-H cleavage TSs ( 3 TS1, 3 TS2, 3 TS5, 3 TS9, and 3 TS8) of 290-300 kJ-mol" 1 
to the initial C-C cleavage TSs ( 3 TS13 and 3 TS16) of -200 kJ-mol 1 . The activation strain appears to 
be related to the bond strength of the activated bond and with the percentage-wise extent of bond 
stretching in the TS. Typical strengths and lengths of the C-H and C-C bonds are the following: 
414 (C-H) and 347 kJ-mol" 1 (C-C), and -1.1 (C-H) and -1.5 A (C-C) [52]. Moreover, we recall that 
the percentage-wise extent of bond stretching in the TS for oxidative insertion is 70-80% (C-H) and 
-30% (C-C). Hence, both the bond strength and the percentage- wise bond elongation in the TS 
decreases from C-H to C-C. This correlates nearly with the activation strain AZ^strain, which decreases 
in the same order. It is indicated that the activation strain prefers the C-C oxidative insertion over the 
C-H oxidative insertion. That is to say, the activation strain AEf stm [ n makes in the reverse order as the 



Int. J. Mol. Sci. 2012, 13 



9293 



activation energy Aff, for the C-H and C-C oxidative insertion. Alternatively, the strength of the TS 
interaction increases from —50 kJ-mol 1 (C-C) to -220—240 kJ-mol 1 (C-H). From C-C to C-H bond 
activation, the strengthening of the TS interaction varies among 170-190 kJ-mol 1 , whereas the 
activation strain changes only by -90 kJ-mol 1 . Thereby, the lowering activation energy from C-C to 
C-H bond activation stems mainly from the strengthening of the TS interaction. 



Table 1. Geometry (in A) at the BPW9 1/6-3 11 ++G(d, p), Lanl2dz level and activation 
strain analysis of transition state (in kJ-mol" 1 ) at the CCSD(T)//BPW9 1/6-3 ll++G(d, p), 
Lanl2dz level for the C-H and C-C bond cleavage in C3H8 by Pt atom. 



Activated 


TS 


Length in 


Length in 


Stretching 


Stretching 


strain 




AE* 


bond 




CsHg 


TS 


in TS 


in TS (in %) 








C-H 


3 TS1 


1.100 


1.947 


0.847 


77.0 


301.2 


-240.0 


61.2 


C-H 


3 TS2 


1.100 


1.872 


0.772 


70.2 


289.6 


-219.9 


69.7 


C-H 


3 TS5 


1.103 


1.948 


0.845 


76.6 


293.3 


-223.8 


69.5 


C-H 


3 TS8 


1.101 


1.879 


0.778 


70.7 


303.8 


-238.3 


65.5 


C-H 


3 TS9 


1.101 


1.948 


0.847 


76.9 


290.1 


-221.5 


68.6 


C-C 


3 TS13 


1.535 


2.032 


0.497 


32.4 


207.5 


-51.4 


156.1 


C-C 


3 TS16 


1.535 


2.032 


0.497 


32.4 


207.5 


-51.4 


156.1 



Thus, the stabilizing interaction AE^[ nt prefers C-H oxidative insertion, whereas the activation strain 
A^^strain favors C-C oxidative insertion. From C-C to C-H oxidative insertion, the lowering of 
activation barrier is mainly caused by the TS interaction AE^ int becoming more stabilizing. 

4. Conclusions 

The reaction mechanism of the gas-phase Pt atom with C3H8 has been systematically investigated 
on the singlet and triplet potential energy surfaces. Considering the initial interaction of Pt atom with 
C3H8, Pt atom prefers the attack of primary over secondary C-H bonds in propane. The major and 
minor reaction channels lead to the dehydrogenation products PtC3H6 + H2 and the deethanization 
products PtCH 2 + C2H6, respectively, whereas the possibility to form the demethanation products 
PtC2H 4 + CH 4 is so small it can be neglected. Over the 300-1 100 K temperature range, the branching 
ratios for the formation of PtC3H 6 + H 2 and PtCH 2 + C2H6 are calculated to be 97.7%~83.6% and 
2.3-16.4%, respectively. The MERP for the formation of the main products PtC3H 6 + H 2 , involving 
one spin inversion, prefers to start at the triplet state and afterward proceed along the singlet state. The 
optimal C-C bond cleavages are assigned to C-H bond activation as the first step, followed by cleavage 
of a C-C bond. This reactivity mode is complementary for the classical reactivity picture through the 
direct C-C cleavage intermediate. 

Furthermore, both the C-H insertion intermediates [l-HPtC 3 H 7 , c/s-l-HPtC 3 H 7 , HPtCH(CH 3 ) 2 , 
3-HPtC 3 H 7 , c/s-3-HPtC 3 H 7 ] and the C-C insertion intermediates [CH 3 PtC 2 H 5 , c/s-CH 3 PtC 2 H 5 , 
C 2 H4Pt(H)(CH3), and CH4PtC 2 H4] are thermodynamically preferred. However, the C-H insertion 
intermediates are kinetically favored over the C-C insertion intermediates. These results are in 
qualitative agreement with the experimental results, in which the C-H insertion product is 
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experimentally observed and the C-C insertion product is not formed in observable quantity in 
Pt + C 2 H 6 system. 

Unlike the role of the activation strain AEf stm { n , the stabilizing interaction AE^[ nt favors the initial 
C-H oxidative insertion over the initial C-C oxidative insertion. From C-C to C-H oxidative insertion, 
the lowering of activation barrier is mainly caused by the TS interaction Alf i nt becoming more stabilizing. 
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Appendix 

Zero-point energies (ZPE) (hartree), total energies (E c ) (hartree) corrected by ZPE, relative energies 
(E r ) (kJ-mol 1 ) of various species at BPW91/6-311++G(d, p), Lanl2dz and 
CCSD(T)//BPW9 1/6-31 1++G(d, p), Lanl2dz levels with respect to the ground reactants Pt( 3 D) + C 3 H 8 . 
The standard orientations and vibrational frequencies of various species calculated at the 
BPW91/6-311++G(d, p), Lanl2dz level in the reactions of Pt + C3H8. The schematic energy diagrams 
in the activation of C 3 H 8 by Pt atom calculated at the CCSD(T)//BPW91/6-311++G(d, p), Lanl2dz 
level, in which the relative energies (kJ-mol 1 ) for the corresponding species relative to Pt( 3 D) + C3H8 
are shown. 
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